Quantum Walks on Sierpinski Gaskets 

Pedro Carlos S. Lara* Renato Portugal} and Stefan Boettcher-'' 
September 11, 2012 



Abstract 

We analyze discrete-time quantum walks on Sierpinski gaskets using 
a flip-flop shift operator with the Grover coin. We obtain the scaling of 
■ two important physical quantities: the mean-square displacement and the 

O | mixing time as function of the number of points. The Sierpinski gasket is a 

fractal that lacks translational invariance and the results differ from those 
described in the literature for ordinary lattices. We find that the displace- 
ment varies with the initial location. Averaged over all initial locations, 
our simulation obtain an exponent very similar to classical diffusion. 

1 Introduction 

Discrete-time quantum walks have been introduced by Aharanov, Davidovich, 
and Zagury [I] as the quantum version of classical random walks. Quantum walks 
on lattices can spread out ballistically, in contrast with the diffusive behavior of 
classical random walks. This characteristic has motivated many studies pursuing 
quantum algorithms that are faster than their classical counterparts [31 0]. 

Discrete-time quantum walks have been investigated previously on many 
graphs. The most studied graph is the one-dimensional line[5l [6l [7]. Quantum 
walks have been analyzed on two-dimensional square lattices [8, 9:, and on the 
hypercube[lO . A spatial search using the discrete-time quantum walk model 
has been undertaken on the Sierpinski gasket [TT]. and on the Hanoi network of 
degree 3. [12] A quantum walk on the dual Sierpinski gasket using the continuous- 
time quantum walk model has been analyzed by Agliari et. al. [13] 

In this paper we focus our attention on discrete-time quantum walks on the 
Sierpinski gasket. We analyze the dynamics based on the standard evolution 
operator U = S ■ (C ® I), where S is the flip-flop shift operator, C is the coin, 
and I is the Identity operator. Throughout, we are using the Grover coin. The 
main physical quantities that we analyze are the mean-square displacement in 
form of the standard deviation in position, the limiting probability distribution, 
and the mixing time. The results are compared with classical random walks 
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on the Sierpinski gasket and with quantum walks on other graphs, such as the 
square lattice. 

The paper is organized as follows: In Section [2] we derive the evolution 
equation for quantum walks on Sierpinski gaskets. In Section [3] we present the 
numerical results for the standard deviation, the limiting probability distribu- 
tion, and the mixing time. In the last section, we present our conclusions. 

2 Standard Quantum Walk Dynamics 

The Sierpinski gasket of generation g is a degree-4 regular graph, an example 
of which is depicted in Fig. [I] for g = 2. It has N = 3(3 9 + l)/2 nodes and a 
fractal or Hausdorff dimension of df = log 3/ log 2, which is larger than for the 
line and smaller than for the plane. 
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Figure 1: Sierpinski gasket of generation g = 2 embedded in a two-dimensional 
plane. 

A coined quantum walk on the Sierpinski gasket embedded in the two- 
dimensional plane has a Hilbert space He 8^?, where He is the 4-dimensional 
coin subspace and Hp is the ./V-dimensional position subspace. Hp is spanned 
by vectors of type \x,y) with integers < x < 2 9+1 and < y < 2 9 restricted 
to be on the gasket, as shown in Fig. [TJ Due to the embedding, we use the 
computational basis {|fc) ,0 < k < 5} for the coin space He, but only four of 
these basis vectors are utilized for each vertex. 

The shift operator for the internal vertices is 

S\k)\x,y) = \-k)\x + f(k),y + g(k)), (1) 

where — k is the inverse of k modulo 6. The coin value is inverted after the shift 
(flip-flop shift). Functions / and g are defined in the Table [21 

For the external vertices (0,0), (2 s , 2 s ), (2 9+1 ,0), where g is the genera- 
tion level, we consider two cases of boundary conditions: (1) periodic and (2) 
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Table 1: Auxiliary functions for the shift operator. 



reflective. In case (1), the action of the shift operator is given by 



S |3> |0,0> 
S |4} |0,0} 
5|1) |2 9 ,2 9 } 
S\2) \2 9 ,2 9 ) 
S\0) \2 9+1 ,0) 
S|5} |2 S+1 ,0) 



-3} |2 9+1 ,0) 
-4} |2 9 ,2 9 } 
-1)10,0} 
-2} |2 9+1 ,0) 
-o) |0,0) 
-5)|2»,2»). 



Those special cases can be implemented through functions / and g and, in this 
case, / and g will depend on the location (x,y). In case (2), the action of the 
shift operator is given by 



S |3} |0,0} 
S |4} |0,0} 
\2 9 ,2 9 ) 
S\2) \2 9 ,2 9 ) 
S\0) \2 9+ \0) 



14} |1,1} 
|3} |2,0} 

|1) \2 9 - 1,2 s - 1} 
|2}|2 S + 1,2 S -1} 
|5) |2 S+1 - 1,1) 
|0}|2 s+1 -2,0). 



S|5)|2» +1 ,0) = 
The Grover coin is defined as 

G = 2|D} (D|-7, 
where |D) = \ J2k=o \^)- ma trix representation is 



G =2 



-1 
1 
1 
1 



1 1 

1 1 

-1 1 

1 -1 



(2) 



(3) 



The generic state of the walker at time t is described by 

5 2 9+1 min{x,2 9 + 1 -rr} 

!*(*)> =EE E Tpk; X ,y(t)\k}\x,y) 

k=0 x=0 y=0 



(4) 



where the coefficients ipk-,x,y(t) are complex functions that obey the normaliza- 
tion condition 

EEi^*.»(*)i 2 = 1 ' ( 5 ) 

fc=0 x,y 
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for all time t. 

Applying the evolution operator 

U = S (G ® 7) (6) 

to the generic state, we obtain 

l*(* + l)> = E E^';^W4? fe ' } \~k)\x + f(k),y + g(k)). (7) 
fc,fc'=0 

Renaming the dummy indices, we obtain 

i*(*+i)> = E E G -Xt ( "* ),, '" fl( " fc)) ^--/(-fc).«-»(-«(*)i fe >i a; .^8) 

Expanding the left hand side of the above equation in the computational basis 
and equating like coefficients, we obtain the evolution equation for the quantum 
walk, 

+ = EGi{, ( - fc) ^ 9( - fc)) ^ ; ,-/(-,)^- s (- fc )W- (9) 

fe'=0 

The matrix G^ x,v ' depends on x, y, since there are six types of vertices that are 
distinct in their orientation. For each one, we have to use the correct labels for 
their edges. 

We use Eq. ([9]) to numerically simulate the evolution of the quantum walk 
using initial conditions of the form |D) \x,y), where |D) = | J2k=o 1^) * s tne 
uniform vector in the coin space. Note that |D) is not biased. The same is 
true for the Grover coin G. This coin and the flip-flop shift operator play an 
important role in spatial search algorithms |17[ ITTj. 

3 Physical Quantities 

In this section, we analyze the behavior of quantum walks on the Sierpinski gas- 
ket with the focus on the diffusion processes. The main physical quantities that 
wc analyze are the mean-square displacement in form of the standard deviation 
in position, and the mixing time. 

3.1 Standard Deviation 

The physical quantities that we will analyze are defined using the probability 
distribution over the vertices of the graph. It is one of the main physical quan- 
tities that is available in the analysis of the behavior of quantum walks. The 
probability distribution is given by 

5 

p(t;x,y) = ElV'fc;^W| 2 - (10) 

k=0 
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The position standard deviation a(t) is defined as 



a{tf 



<J x (t) 2 + a y {t) 2 , 



(11) 



where 



Ox{t) 2 = ^x 2 p(t,x)-[^xp{t,x) 



X \ X 



v \ y 




and 



p(t,x) 



y 



p(t,y) 
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For a walker that starts located on a specific vertex, the standard deviation 
at intermediate times 1 <C t <C t co increases as a power-law 



which defines the diffusion exponent 1 < d w < oo, in analogy to a classical 
walk. On a finite system, the walk eventually reaches the farthest vertex which 
cuts off the growth in a at some time t co , beyond that it oscillates around an 
average value. In order to analyze the diffusion process, we are interested in the 
asymptotic behavior of the power law regime for the infinite system, g — > oo. 
Therefore, on any finite Sierpinski gasket, we bound the evolution time to be 
smaller than the time that the walker takes to reach the farthest vertex. To mea- 
sure displacement (standard deviation), we use reflective boundary conditions 
and the flip-flip shift operator. 

The first result we have obtained from the simulations, which is strikingly 
different from the behavior of quantum walks on lattices, is that the scale of the 
standard deviation depends on the initial vertex. For example, for generation 
g = 8, the fastest growth in the displacement is obtained when the walker starts 
on vertex (x — 247, y = 5) Fig. [5] shows the standard deviation as function 
of the number of steps of a quantum walk with the initial state |D) 1 247, 5) . 
Our fit yields a = Lit 52 , i.e., d w s» 1.92, significantly faster than classical 
diffusion (d w = log 2 5 = 2.32 . . .) but still considerably slower than the standard 
deviation for quantum walks on a square lattice, for which d w = 1. 

The slowest-growing standard deviation of a flip-flop quantum walk on the 
Sierpinski gasket of generation g = 8 is obtained for the initial state |D) 1 256, 224) . 
The best fit in this case using the data in Fig. [2] is a ~ 2.0i 29 (or d w w 3.45), 
which shows a very small spreading rate characterizing a sub-diffusive process. 



a(t) 



(12) 
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Figure 2: Standard deviation rr(t) of a flip-flop quantum walk on the Sierpinski gasket 
(g = 8) after 128 steps: (1) with initial state |D) 247, 5) displayed with circles, (2) 
with initial state |D) 256, 224} displayed with triangles, and (3) averaged over all initial 
states |D) \x,y), \/(x,y) displayed as squares. 



Since the displacement apparently depends on the initial vertex, it is in- 
teresting to define a mean standard deviation a(t) as function of time in the 
following way 

*(*) = ^ !>».»(*)> ( 13 ) 

x,y 

where sub-indices x, y of cr x ,y(t) indicate the initial location used to obtain the 
standard deviation. From each location the walker evolves from an initially 
uniform state in coin space. Fig. [5] depicts the behavior of <r(t). The numerical 
results suggest a best fit of a(t) ~ 1.3t 44 , or d w ~ 2.27. Average over all 
initial locations makes this scaling exponent d w a characteristic of the Sierpinski 
gasket of generation g = 8, which happens to be remarkably close to the result 
for classical diffusion, d w = log 2 5 = 2.32 . . .. The histogram in Fig. |3] shows 
the number of such initial conditions that have the same fitted exponent, d w , in 
Eq. (|12p . For example, for d w = 2.14 there are around 450 vertices that can be 
used as initial condition to obtain the same scaling. The range in scaling (for 
g = 8) extends from d w ~ 1.92 to 3.45, although the bulk of the distribution is 
centered very close to the result for classical diffusion, marked by a vertical line. 
The width of this distribution is only about 17% of the mean, and it would be 
interesting to see whether the width narrows further for increasing system sizes 
g -> oo. 
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Figure 3: Histogram of the frequency of the initial conditions that have the same scale 
in the standard deviation of a flip-flop quantum walk on the Sierpinski gasket at g = 8 
with the initial state |D) \x,y), \/(x,y). A vertical line indicates the corresponding result 
for classical diffusion. 

3.2 Limiting Distribution 

The average probability distribution is given by 

1 T_1 

p(T,x,y) = — 22p(t,x,y). (14) 
t=o 

Note that p(T, x, y) is a probability distribution for all T, because 

N-l 

p( T > x >y) = !■ 

x,y— 

The interpretation of p(T, x, y) uses projective measurements, therefore p(T, x, y) 
evolves stochastically, and converges to a limiting distribution when T goes to 
infinity. The definition of the limiting probability distribution is 

n( x > v) = J im p( t > x i v)- ( 15 ) 

This limit exists and can be calculated explicitly if the expressions for the eigen- 
values of the evolution operator are known. The limiting distribution depends 
on the initial condition in general. 

Fig. [4] shows the limiting distribution as a function of position a; of a flip-flop 
quantum walk that departs from the central bottom vertex (x — 256, y = 0) with 
periodic boundary conditions. The probabilities in y-direction have been added 
up to generate a one-dimensional plot. In general, the numerical simulations 
show that the limiting distribution for the Sierpinski gasket depends on the 
initial condition and is highly concentrated around the initial vertex. Note 
that a walker encounters frequent bottlenecks that inhibit spreading. If the 
walker starts in vertex (256,0), there are only two passage points, (128,128) 
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Figure 4: Limiting distribution of a flip-flop quantum walk on the Sierpinski gasket 
(g = 6) with the initial state |D) 64, 0} as function of position x. We have added 
all probabilities with different values of y having the same value of x. The boundary 
conditions are periodic. 

and (384, 128), towards the top of the Sierpinski gasket (in Fig. [T] these points 
correspond to (2, 2) and (6, 2)). A three-dimensional plot shows that the limiting 
probability is very close to zero when y > 1. 
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Figure 5: Total variation distance between p(T, x, y) and ir(x, y) for a flip-flop quantum 
walk on the Sierpinski gasket (g = 6) with the initial state |D) |64, 0). The best fit is 
||p(T, x, y) — tt(x, y)\\ = 1.76/T. The inset show the plot without using the log-log scale. 

The average distribution p(T, x,y) converges to the limiting distribution 
Tr(x,y). This can be confirmed by the graph of the distance between these 
distributions as function of time. The total variation distance between two 
probability distributions p{x,y) and q(x,y) is defined as 

II p ~ q H = \p( x >y) ~q(x,y)\ ■ (16) 
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Fig. [5]shows the plot of \\p(T, x, y) — tt(x, y)\\ as function of the number of steps 
represented by T. The best fit suggests that this distance scales as 1/T. We note 
that it decays approximately as 1 /T not only for the Sierpinski gasket but also 
for two-dimensional lattices[5] and hypercubes|10|. This numerical result can 
be explained analytically. Expanding the quantum walk state in the eigenbasis 
that diagonalizes the evolution operator, one can obtain an explicit expression 
for \\p(T,x,y) — n(x,y)\\ for any initial condition. The expression has the form 

\\p(T, x, y) - n(x, y)\\ = ±2& ^ H, (17) 

where A c is a difference between two non-equal eigenvalues, c x . y is a constant, 
and the inner sum is over all pairs of non-equal eigenvalues. Some results of 
Aharonov et. al.[T5] help to obtain that analytical expression. The modulus 
of the term e 2, ™ AcT — 1 in the numerator of \\p(T, x, y) — ir(x, y)\\ is a bounded 
oscillatory function. So, the distance between the average and the limiting distri- 
bution scales as 1/T in general and oscillates around the curve 1/T, confirming 
the data shown in Fig. [5] 



3.3 Mixing Time 

The mixing time r £ is defined as 

r e =min{T|Vt > T, || p(t, x, y) - w{x, y) ||< e}, (18) 

which can be interpreted as the smallest number of steps such that the distance 
between the average distribution and the limiting distribution becomes perma- 
nently smaller than e. If \\p(T,x,y) — n(x,y)\\ obeys an inverse power law as a 
function of time, then r e obeys an inverse power law as a function of e. 

The mixing time depends on the initial condition in general and on the size 
N of the graph. We have generated the same kind of data of Fig.[5]for Sierpinski 
gasket of generation 6 up to 10. Using the best-fitting curves we can estimate 
r e as a function of N. Fig. |6] shows that r e has a power law in terms of the 
number of vertices. The data allows us to estimate that 

r c = i^—j (19) 
when we take |D) |2 9 ,0) as initial state. 



4 Conclusions 

We have analyzed the flip-flop discrete-time quantum walk on the Sierpinski 
gasket of finite generation embedded in the two-dimensional plane using reflec- 
tive and periodic boundary conditions. Our investigations focus on the following 
physical quantities: the position standard deviation (with reflective boundary 
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Figure 6: The mixing time of a flip-flop quantum walk on the Sierpinski gasket with 
the initial state |D) 2^,0). The best fit is r e = 0.034A r0 - 54 . The first point corresponds 
to g = 6. 

conditions) and the mixing time (with periodic boundary conditions) . Perform- 
ing numerical simulations on Sierpinski gaskets up to generation g — 10, we 
have obtained the scaling exponent of the standard deviation as function of the 
number of steps and the scale of the mixing time as function of the number of 
vertices. 

For the system sizes studied, the results depend significantly on the initial 
condition. As fractal lattices lack translational invariance, quantum interference 
effects likely vary strongly with the initial location. A characteristic way to 
assign a distinct diffusion exponent d w to the Sierpinski gasket is provided by 
performing an average over all initial locations. In this case, we obtain an 
average exponent d w ps 2.27 that is remarkably close to the result for classical 
diffusion, d w = log 2 5 = 2.32 . . ., on this system. Therefore, a quantum walk on 
the Sierpinski gasket spreads slower than on a square lattices. 

The limiting distribution for the Sierpinski gasket depends on the initial 
condition and is concentrated around the initial vertex. This happens with 
other graphs such as the two-dimensional lattice 0. The scaling of the mixing 
time for the Sierpinski gasket O(N M /e) is close to the scaling of mixing time 
for the two-dimensional lattice, which is believed[9] to be O^N log N /e). Our 
data is not precise enough to determine the presence of a term that depends on 
logiV. The result differs from the scaling on the cycle which is believed [18] to 
be 0{N\ogN/e). 
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